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We investigate a calculation method for the gravitational evolution of an extreme mass 
ratio binary, i.e. a binary constituting of a galactic black hole and a stellar mass compact 
object. The inspiralling stage of this system is considered to be a possible source of de- 
tectable gravitational waves. Because of the extreme mass ratio, one may approximate such 
a system by a black hole geometry (a Kerr black hole) plus a linear metric perturbation 
induced by a point particle. With this approximation, a self-force calculation was proposed 
for a practical calculation of the orbital evolution, including the effeqt of the gravitational 
radiation reaction, which is now known as the MiSaTaQuWa self- force EP . In addition, a radi- 
ation reaction formula was proposedEf as an extension of the well-known balance formula of 
Press and TeukolskyEP. The radiation reaction formula provides us a convenient method to 
calculate an infinite time averaged loss of "the constants of motion" (i.e. the orbital energy, 
the z component of the angular momentum, and the Carter constant) through the gravita- 
tional radiation reaction, with which one may approximately calculate the orbital evolution. 
Because these methods are approximately equivalent, we investigate the consequence of the 
orbital evolution using the radiation reaction formula. 

To this time, we have used the so-called adiabatic approximation of the orbital evolution 
and considered a method to evaluate the MiSaTaQuWa self-force by use of a linear metric 
perturbation. In this approach, we point out that there is a theoretical question concerning 
the choice of the gauge condition in the calculation of the MiSaTaQuWa self-force. Because 
of this gauge ambiguity, there is a case in which the MiSaTaQuWa self-force might not pre- 
dict the orbital evolution in a physically expected manner, and this forces us to calculate 
a waveform only in the so-called dephasing time. We discuss the reason that such an un- 
expected thing happens and find that it is primarily because we consider the linear metric 
perturbation separately from the orbital evolution due to the self-force. 

We propose a new metric perturbation scheme under a possible constraint of the gauge 
conditions in which we obtain a physically expected prediction of the orbital evolution caused 
by the MiSaTaQuWa self-force. In this new scheme of a metric perturbation, an adiabatic 
approximation is applied to both the metric perturbation and the orbit. As a result, we are 
able to predict the gravitational evolution of the system in the so-called radiation reaction 
time scale, which is longer than the dephasing time scale. However, for gravitational wave 
detection by LISA, this may still be insufficient. We further consider a gauge transformation 
in this new metric perturbation scheme, and find a special gauge condition with which we 
can calculate the gravitational waveform of a time scale long enough for gravitational wave 
detection by LISA. 



§1. Introduction 



Constructing a method for the accurate calculation of gravitational waveforms 
is an essential step toward success in gravitational wave detection. Since the Ein- 
stein equation is nonlinear, it is difficult to derive fully analytic solutions of every 
gravitational wave source, and we therefore consider an approximation method for 
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a specific astrophysical source. 

In this paper, we consider an extreme mass ratio binary system. X-ray and a 
variety of other observations suggest that there may exist supermassive black holes 

(M ~ 10 6 10 s Mq) at the centers of galaxies. A number of stars have been 

observed at the galactic centers, and they interact with each other gravitationally. 
As a result of this multi-body interaction, we may expect that some galactic black 

holes may capture a stellar mass compact object (m ~ 1 lOOM©) deeply into its 

gravitational potential without tidal disruption. We refer to such a pair of objects as 
an extreme mass ratio binary system, and such systems are believed to be possible 
sources of gravitational waves detectable by LISA&. Due to its orbital energy and 
angular momentum, the stellar mass compact object inspirals around the black hole 
for a long time, referred to as the 'inspiral stage' of the binary evolution. During 
this stage, the system has a huge quadrupole moment and emits detectable gravita- 
tional waves, while the system slowly loses its orbital energy and angular momentum 
through this emission of gravitational waves. 

We consider a method to compute the gravitational waveform from the inspiral 
stage for such an extreme mass ratio binary. Because the black hole dominates the 
gravitational potential due to its huge mass, we assume that a metric perturbation 
provides a good approximation to describe the gravitational evolution of the system. 
We employ a Kerr black hole as the background because of the black hole uniqueness 
theorem in general relativity. Since the size of the stellar mass compact object is 
much smaller than the background curvature, it can be treated as a point particle 
to induce a metric perturbation U>. In this picture, the evolution of the binary can 
be interpreted as the orbital evolution of the particle around the black hole. 

It has been asserted that, because we already have a convenient formula to 
calculate a linear metric perturbation in a Kerr background a\ the remaining issue 
in this problem is to develop a method to compute the orbit of a point particle. The 
orbital evolution can be derived by calculating the self-force, and a general formula 
for the self-force employing a linear metric perturbation scheme has been derivedtP . 
Because the linear metric perturbation induced by a point particle diverges along 
the orbit, the implementation for the self- force involves a regularization calculation. 
We have made good progress in obtaining an exact derivation of the self- force & by 
implementing this regularization calculation. However, there remains a problem with 
this method when we consider the self-force of a general orbit in a Kerr background. 
A related method that employs the so-called radiation reaction formula has also been 
proposed a. With this method, one can approximately derive the self- force, although 
it only contains the radiative part of the self-force. However, this method does not 
involve a complicated regularization calculation, and the numerical method has been 
established even for a general orbit in a Kerr background. 

Although the calculation of the self-force seems to be a reasonable step for a 
gravitational waveform calculation, the self-force is entirely gauge dependent and 
there is a question of what physical content the self- force actually has. In fact, one 
can easily show that, with an extreme choice of the gauge condition, the self-force can 
be made to vanish over the entire time interval for which a usual metric perturbation 
scheme is valid in describing a gravitational evolution of the system. We discuss that, 
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although this is mathematically allowed in the usual metric perturbation scheme, 
the resulting description is technically disadvantageous from an energetics point of 
view. If we adopt such a gauge condition, we may eventually have to consider the 
gravitational evolution of the system in a nonperturbative manner, and we must use 
numerical relativity to calculate the nonlinear metric perturbation, to our knowledge. 

Instead of numerical relativity, we consider the extension of the usual metric 
perturbation with a certain constraint on a gauge condition. A primary reason for 
such an unexpected consequence in the usual metric perturbation scheme is that the 
particle must move along a geodesic to the leading order of the perturbation. As a 
result, a linear metric perturbation cannot track gravitational evolution together with 
the orbital evolution affected by the self-force. We propose a new metric perturbation 
scheme in which the MiSaTaQuWa self-force does not vanish. In this approach, 
the orbital evolution caused by this self-force turns to be essential to calculate the 
gravitational waveform. 

The new metric perturbation scheme may be able to predict the gravitational 
evolution of the system over a time scale longer than that of the usual metric per- 
turbation. However, this time scale may not be sufficient to calculate a gravitational 
waveform for LISA. For this reason, we further consider a gauge transformation along 
an orbit so as to extend the validity of this new metric perturbation scheme. We 
find a special gauge condition with which the self-force contains only the radiation 
reaction component of the system, and, in this gauge condition, the adiabatic ap- 
proximation of the linear metric perturbation can predict the evolution over a time 
scale long enough for LISA. 

This paper is organized as follows. In Sec.||, we consider the gauge ambiguity of 
the MiSaTaQuWa self-force in the usual metric perturbation scheme and clarify the 
reason why such a prolblem could arise in the nonlinear theory of radiation reaction. 
In Sec|| we propose a new metric perturbation scheme, with which one can treat the 
gravitational evolution of the system together with the orbital evolution. In Sec.[|, we 
consider a further gauge transformation in the new metric perturbation scheme and 
propose the so-called radiation reaction gauge condition. We summarize our result 
in Sec.|5|, with emphasis on clarifying the validity of this new metric perturbation 
scheme. 

§2. Gauge ambiguity in the self- force problem 

Gravitational radiation reaction is a physically real phenomenon, since we can 
define a gauge invariant momentum flux of gravitational waves in an asymptotic flat 
region of a background metric However, we find that a self- force is entirely gauge 
dependent consistently with a usual metric perturbation scheme and that it can be 
made to vanish through a special gauge choice along an orbit over the entire time 
scale for which the usual metric perturbation scheme is valid. 

It is well-known that the self-force can be made to vanish at an orbital point 
by making the appropriate gauge choice, but we conjectured this might not be a 
problem, because the self-force would have non-vanishing components if we take 
the average of the self-force over a long time. In Ref. we show that this is true 
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if we take the average of the self-force over an infinite time-J, and we prove that 
non-vanishing components are gauge invariant in agreement with the gravitational 
radiation reaction. However, in the usual metric perturbation scheme, one cannot 
consider an infinite time average of the self-force because it is valid only over a finite 
time interval. 

In Subsec.2.1, we consider a self-force in the usual metric perturbation scheme, 
and, in this scheme, we find that the average of the self-force can be made to vanish 
by making a special gauge choice along the orbit. This suggests that the relation 
between the self-force and the gravitational radiation reaction is entirely gauge de- 
pendent, and that it is not trivial as we previously believed. As a result, the so-called 
adiabatic approximation of the orbitLP'tf would yield a gauge dependent prediction 
for the orbit. Because we expect that the orbital evolution caused by the self-force 
can be observed by the modulation of gravitational waves, a physically reasonable 
prediction of the orbit must be gauge invariant. This leads to the questions of 
how the self-force is related to the gravitational radiation reaction and what gauge 
conditions we should use to calculate the self-force. We discuss these questions in 
Subsec.2.2 

2.1. Metric perturbation 

We assume that a regularization calculation for the self-force is formally possible 
by using a matched asymptotic expansion of the metric of the near-zone expansion 
and the metric of the far-zone expansion, as considered in Ref. Because the 
resulting self-force is defined by regularizing the metric of the far-zone expansion, 
we consider only the metric of the far-zone expansion. The metric of the far-zone 
expansion consists of the sum of a regular vacuum background metric and its metric 
perturbation induced by a point particle moving in the background metric. We 
assume the point particle has an appropriate internal structure™ so that the metric 
of the far-zone expansion matches the metric of the near-zone expansion. 

Suppose we use the usual metric perturbation scheme for the calculation of 
the metric of the far-zone expansion. Then the far-zone expansion is defined as a 
perturbation expansion with the small parameter m/L where m is the Schwarzschild 
radius of the particle and L is the curvature scale of the background along the orbit. 
In the usual metric perturbation scheme, we expand the metric and the stress-energy 
tensor in term of this small parameter as 

9,u = g% 9) + (m/L)h$ + {m/Lfh$ + ■ ■ ■ , (2-1) 
T" v = {m/L)T^ v + (m/L) 2 T^ u + • • • , (2-2) 

where gffl is the vacuum background metric. For a valid perturbation, we must 
have 

O(l) > \{m/L)h$\ > \{m/L) 2 h$\ > ■•• , (2-3) 
O(l) > \{m/L)T^ v \ > \{m/L) 2 T^ u \ >■■■ . (2-4) 

*' It is important to note that this argument depends only on the properties of the orbital 
equation, and it does not depend on the metric perturbation scheme. 
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We also expand the Einstein equation in m/L, and schematically we have 

£(1)^(1)] = r (i) ( * ) (2-5) 

^(2)] + [ h w > = t (2)f ; ( 2 . 6 ) 



where and G^ 2 ^ v [h, h] are respectively the terms linear and quadratic in 

h^v of the Einstein tensor G^ v [g + h] . 

In Ref.B*, it is found that T ( - 1 ^ u is the usual stress-energy tensor of a monopole 
particle for a matched asymptotic expansion of the metric. It is important to note 
that the linearized Einstein equation is algebraically divergence free with respect to 
the background metric, i.e. G^ ,xv - v = T^^ u - U = 0. Thus, in this metric pertur- 
bation scheme, the particle must move along a geodesic in the background metric 
gffi for a consistent solution of (2-5). The explicit form of T^ 2 ^ u can be derived 
by carrying out a further matched asymptotic expansion of the metric, but this 
has not been obtained yet. However, at least, it is clear that it has a term of the 
monopole particle coming from the deviation of the orbit from a geodesic because 
the MiSaTaQuWa self-force can be derived from LHS of (2-6) through a mass renor- 
malizationi^. This shows that in the expansion of the stress-energy tensor (2-2), we 
assume an expansion of the orbit by m /L as 

zHt) = z^(t) + (m/L)z^ + ■■■ , (2-7) 

where r is defined to be the proper time in the background metric, and z {b s^ is a 
geodesic of the background metric given by T^'^ u - V = 0. For a valid perturbation, 
we must have 

0(1) > \(m/L)z^\ > ■■■ . (2-8) 

Intuitively, the orbit deviates from a geodesic as a result of the gravitational radiation 
reaction, and eventually the condition (2-8) will be violated. Then, beyond this time, 
this perturbation scheme would fail to approximate the system. 

We consider a gauge transformation over a time interval during which the metric 
perturbation still approximates the system. The gauge transformation is defined 
as the small coordinate transformation — > = + (m/L)^(x). The orbit 
transforms as z^(t) — > z^(t) = z^(t + (m/L)Jr) + (m/L)^(z(r)), where 5t is 
a function of r, so that r remains the proper time of the new orbit z^{t) in the 
background metric. Applying this to the perturbation expansion of the orbit (2-7), 
we obtain the gauge transformation of (m/L)z^^ as 

(m/L)z^ l (T) -> (m/L)z^( T ) 

= (m/L)z^(T) + (m/Ljv^Sr + (m/L)^(z^ (t)$-9) 

where v^ b9 ^ = dz^ b9 > /dr. We can eliminate (m/L)z ( - 1 ^(T) if the gauge transforma- 
tion along the orbit satisfies 



{m/L)e{z (bg H-r)) = -{m/L)z {1)fl {T) + v^{m/L)A T , (2-10) 
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with an arbitrary function At. There always exists a gauge transformation which 
satisfies this condition over the entire time interval during which the metric pertur- 
bation scheme is valid, because it is smaller than O(l) in accordance with (2-8), and 
the orbit becomes a geodesic of the background metric under this gauge condition. 

Because the MiSaTaQuWa self-force causes a self-acceleration of the leading or- 
der in m/L-expansion and causes the orbit to deviate from the geodesic, the relation 
(m/L)z^ 1 ^(r) = implies that the MiSaTaQuWa self-force entirely vanishes under 
this gauge condition over the entire time interval for which this metric perturbation 
scheme is valid. This extreme example suggests that the MiSaTaQuWa self-force 
is completely gauge dependent in this usual metric perturbation scheme, and, even 
the time average of the self-force over a long time scale could be gauge dependent in 
general, as long as we use this perturbation scheme. 

At this stage, we see the mathematical reason why the self-force may not include 
the effect of gravitational radiation reaction. Because this perturbation scheme only 
allows a small deviation from the background geodesic, as expressed by (2-8), one 
could always eliminate this deviation by a gauge transformation, as in (2-9). We note 
that this problem cannot be solved by calculating a non-linear metric perturbation 
because the problem resides in the metric perturbation scheme itself. This shows 
that we must modify the metric perturbation scheme so that we can describe a non- 
perturbative orbital deviation from the background geodesic as (m/L)^^^) ~ 
0(1). If the orbit can deviate from the geodesic nonperturbatively, it cannot be 
eliminated by the gauge transformation. 

This extreme example suggests that the relation between the self-force and the 
gravitational radiation reaction is not trivial, as we believed before&'H*. By an 
adiabatic approximation of an orbit, we usually use the time averaged part of the 
self-force. If this is a component averaged over an infinite time, it is known to be 
gauge invariant, and we have a gauge invariant prediction for the orbit. However, 
because the metric perturbation is not valid over an infinite time scale, it is not valid 
to consider an infinite time average, and it seems reasonable to take a finite time 
average. Then the orbital prediction becomes gauge dependent. In this situation, 
if we still use a linear metric perturbation, the gravitational waveform at infinity 
would also be gauge dependent, although we would expect it to be gauge invariant 
because it is observable. 

In the next subsection, we study this example from a different point of view, 
using a theory of a nonlinear metric perturbation. We expect that this extreme 
example may reveal to us the nature of the self-force, and, consequently we will be 
able to understand why and how we should take a certain class of gauge conditions 
for the self-force. 

2.2. Energetics of the orbit and radiation 

Some idea of why such an unexpected thing happens can be gained from Ref.ik 
In Sec. Ill of Ref.@) , we argue that the orbital energy does not decrease monotonically 
as a result of the emission of gravitational waves. We conjecture that gravitational 
radiation has its own energy and that the self-force describes the interaction between 
the orbital energy and the radiation energy. In fact, one can define an effective stress- 
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energy tensor for the orbit and gravitational radiation, as argued in Ref.Q). As in the 
previous subsection, we consider only the metric of the far-zone expansion, and we 
suppose that the metric consists of the sum of the vacuum background metric and 
its perturbation induced by a point particle with an appropriate internal structure EP 

The Einstein equation for the far-zone expansion can be written formally as 
G^}[h] + G^\h] = Tfj, v , where G^} and G^t + ^ contain the linear terms and the rest 
of the Einstein tensor with respect to h^ u , respectively. Because G$ is algebraically 
divergence free, one can define a conserved stress-energy tensor in the background 
metric as 

T>™ = - G (2+) H/i] . (2-H) 

We consider that the first and second terms of (2-11) represent the effective stress- 
energy tensors for the orbit and the gravitational radiation, respectively, and neither 
of these is conserved in the background metric by itself. Although the effective stress- 
energy tensor (2-11) is defined in a nonlinear manner in a general curved background, 
it reduces to the well-known gauge invariant stress-energy tensor for gravitational 
wavesin the so-called short wavelength approximation and in the weak perturbation 
limit B> . 

We suppose that the background metric g^S is a Kerr black hole and that it 
has a Killing vector £„. Then we can define total, orbital and radiation energies in 
the background metric as 

E (tot) _ fi(orb) _|_ E (rad) 

E (orb) = j d £^ vT W ^ E (rad) = _ J dS ^ G (2+)^ > ( 2 -12) 

where the surface integral is taken over the spacelike hypersurface bounded by the 
future horizon and the future null infinity of the background black hole metric, and 
we assume an appropriate regularization along the orbit of the particle, say, by a 
matched asymptotic expansion. We note that the radiation energy 

E (rad) 

is not 

necessarily positive. 

Integrating (2T1) over a smalL world tube surface around the orbit, one can 
derive the MiSaTaQuWa self- forced. This shows that the self-force describes the 
interaction between the orbital energy and the gravitational radiation energy. The 
gauge ambiguity of the self-force could be interpreted as an ambiguity in defining the 
orbital energy and the radiation energy of (2-12), since and G( 2 +)^[/i] are gauge 
dependent. In order to see this, we apply the metric perturbation scheme presented 
in Subsec.2.1 to (2-12). This yields 

= ( m /L)E(° r W + (m/L)E {orb ^ + ■■■ , 
{m/L)E^ orb ^ = J dU^ v {m/L)T^ u , • • • (2-13) 

E (rad) = ( m / jL )2 jB (rad)(2) + ... ^ 

(m/L) 2 E {rad ^ = - J dU^(m/L) 2 G {2)flu [h^\h^} , • • • . (2-14) 
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In the usual perturbation scheme, (m/L)T^ 1 ^ u becomes the conserved stress-energy 
tensor of a monopole particle moving along a background geodesic, and £ , ( orfc )( 1 ) 
becomes a constant of motion. (m/L) 2 T^ fll/ has not yet been derived explicitly, 
and calculations of (m / L) 2 i£( or6 )( 2 ) and {m/ L) 2 E^ rad ^ 2 ^ would involve regulariza- 
tion provided by the matched asymptotic expansion. However, it includes a contri- 
bution from the orbital deviation of the monopole particle relative to the background 
geodesic described by (m/L)T^ tJ,u , as we discussed in the previous section. A gauge 
transformation of this part is defined by £( m / L)^( m I 'V)T^ lUJ , and we find that the 
orbital energy transforms under gauge transformations as 

= (m/L) 2 E (orb W + (m/L) 2 5E , (2-15) 

where 5E is arbitrary. Thus, the orbital energy is gauge dependent in the usual 
metric perturbation scheme. 

By contrast, the balance formula considers the gauge invariant momentum flux 
of gravitational waves at the future null infinity and the future horizon of the back- 
ground black hole geometry. Hence, it describes the radiation reaction to the total 
energy, E^ tot \ rather than that to the orbital energy, 

Although we have not yet proposed a new metric perturbation scheme in which 
a metric perturbation would be valid over a much longer time scale (see the next 
section), the orbital energy and the radiation energy are defined in a nonlinear man- 
ner as (2-12), and we cany consider what happens when we keep the orbital energy 
constant under a gauge transformation while the total energy continues to decrease 
through the gravitational radiation reaction over a sufficiently long time. Because 
the radiation energy, E^ rad ^ , is not necessarily positive, it is easily seen that the am- 
plitude of the metric perturbation would inevitably increase, because the radiation 
energy, E^ rad > , would decrease as a result of the radiation reaction. Although this is 
mathematically allowed, it creates two technical problems: 1) A linear approxima- 
tion can no longer be used in calculating the evolution of the system gravitationally. 
2) The background coordinate system cannot be used as an approximate reference 
to measure the orbit. 

The above considerations reveal what is happening with regard to the behavior 
discussed in the previous subsection. It is reasonable to use a finite time averaged 
self- force with an arbitrary gauge condition for the adiabatic evolution of the orbit, 
and we would then have a gauge dependent prediction for the orbit. However, it 
is not reasonable to use a linear metric perturbation to calculate the gravitational 
waveform at infinity if we do not make an appropriate choice of the gauge condition 
for the self-force. We believe that the gauge invariant nature of the gravitational 
waveform at infinity can be recovered with an appropriate nonlinear calculation for 
the metric in this case. 

*' The radiation reaction in the balance formula is gauge invariant B, and one may expect that 
the total energy is also gauge invariant. However, an explicit computation to prove its gauge invari- 
ance is not trivial because the definition of the total energy involves a gauge dependent regularization 
calculation and a robust volume averaging operation. 
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In order to avoid the difficulty involved in dealing with a nonlinear metric pertur- 
bation, it is disadvantageous to use such a gauge condition, and we therefore assume 
consider that the orbital energy must decrease through the radiation reaction for a 
reasonable choice of gauge conditions, so that one can understand the orbit in the 
background coordinate system as an approximate reference. We define the averaged 
energetic gauge condition along the orbit by 

where <> represents a long-time average. 

Although the above condition is a necessary condition for a reasonable gauge 
choice, this may not be a sufficient condition, because the orbit in a Kerr back- 
ground is not characterized by the energy E^ orb ^ alone. In order to confirm that 
a certain gauge condition is actually advantageous, it is necessary to demonstrate 
that we can actually construct a new metric perturbation scheme that is valid over 
a longer time scale. For this, the new metric perturbation scheme must allow a 
nonperturbative deviation of an orbit from a geodesic, and therefore it would not be 
constrained by (2-8). In the next section, we propose an adiabatic approximation of 
the metric perturbation in which the usual metric perturbation scheme is extended 
adiabatically. 

§3. Adiabatic approximation 

The adiabatic approximation of the orbital evolution is well known and is used 
in Ref. lP'EP. It consists of the following two steps: 

1 ) From the fact that the self-force is weak, we can instantaneously approximate the 
orbit by a geodesic, and, instead of the orbital coordinates, we consider the evolution 
of geodesic "constants" . 

2) We calculate the evolution of the geodesic "constants" caused by the gauge invari- 
ant momentum flux of gravitational waves induced by the approximated geodesic. 

In a Kerr background, Step 2) is insufficient because we do not have a gauge invari- 
ant radiation reaction to the Carter "constant". Instead of Step 2), the radiation 
reaction formula @) includes the following step: 

2 ') We calculate the evolution of geodesic "constants " employing the gauge invariant 
"infinite time av era ged" self-force induced by the approximated geodesic. 
As shown in Ref.li^P, the "infinite time averaged" self- force acting on the orbital en- 
ergy and the angular momentum is identical to the gauge invariant momentum flux 
of gravitational waves. Hence, 2') appears to be a generalization of 2). 2') seems 
more reasonable than 2) because 2') is based on the self-force. However, one cannot 
take an infinite time average in the usual metric perturbation scheme because the 
metric perturbation scheme is valid only in an finite time. Instead of 2'), we can 
consider the following step: 

2") We calculate the evolution of geodesic "constants" employing the self- force in- 
duced by the approximated geodesic. 

In this case, the self-force is entirely gauge dependent, and the orbital evolution 
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becomes gauge dependent as argued in the previous section. 

Because general relativity allows an arbitrary coordinate system, an orbit is 
always "gauge dependent". However, in a metric perturbation scheme, we use a 
fixed coordinate system of the background geometry, and we can define a gauge 
invariant object using the background metric as an approximate reference. To have 
a gauge dependent prediction for the orbit resulting from the self-force means that the 
coordinate system of the background geometry no longer acts as a proper reference 
for the full geometry. As shown in the previous section, this happens if the system 
cannot be described properly by a linear metric perturbation with an inappropriate 
gauge condition. This suggests that any argument for a self-force contains a gauge 
ambiguity until we consider the problem together with the evolution of the metric 
perturbation, so that we can carefully consider the validity of the metric perturbation 
during the orbital evolution. For this purpose, it is necessary to have a new metric 
perturbation scheme which predicts the orbital evolution consistently without the 
orbital constraint (2-8). 

In Subsec.3.1, we propose an adiabatic approximation of the linear metric per- 
turbation. In Subsec.3.2, we evaluate a metric perturbation in this new metric 
perturbation scheme based on the radiation reaction formalism!^. In Subsec.3.3, we 
investigate the orbital evolution and consider the validity of this metric perturba- 
tion. The reason for using the radiation reaction formula is that we have an explicit 
general form of the self-force in the radiation reaction formula, and we can identify 



gauge invariant "infinite 
. (see (3-14)) Thus, one 



the component of the self-force which corresponds to tt 
time averaged" momentum flux of gravitational wavest 
can easily formulate a new metric perturbation scheme under the averaged energetic 
gauge condition (2T6) using this explicit form. 



3.1. Adiabatic linear metric perturbation 

As we summarize in Appendix [A], an inspiralling geodesic in a Kerr geometry is 
conveniently characterized by 7 constants with an orbital parameter A. We denote 
the orbital constants by 7 = {£ a ;a = E, L,C; X b ;b = r,6;C c ;c = t,<fi}. (See Ap- 
pendix [A| for definitions.) In the usual metric perturbation scheme, a linear metric 
perturbation is induced by a particle moving along a geodesic 7. Therefore, it is a 
function of 7, and we write /$$(x;7). 

We consider extending this linear metric perturbation by including the effect of 
the self-force. Due to the self-force, the "geodesic constants" evolve as functions of 
the orbital parameter A. We foliate the background spacetime into smooth spacelike 
hypersurfaces that intersect the orbit M We define the foliation function f(x) such 
that it is equal to the orbital parameter at the intersection with the orbit: f(z(X)) = 



*' The condition that the foliation surfaces are space-like is not important, as long as we consider 
the adiabatic approximation of the linear metric perturbation. For a well-defined nonlinear metric 
perturbation, there must be appropriate attenuation behavior of the adiabatic metric perturbation 
at the past null infinity. Otherwise, an asymptotic regularization is needed to calculate the nonlinear 
metric perturbation. Eif 
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A. Using the foliation function, we define an adiabatic linear metric perturbation as 

h^ 1 Hx)=hW(x-Mf(x))). (3-1) 

Because there is an ambiguity in the foliation function, the adiabatic metric per- 
turbation is not defined uniquely. However, our definition is sufficient for a leading 
correction to the gravitational evolution of the system. 

Because the derivative acts on f(x), we have [h^.x( x ',7)]^=-y(f(x)) 7^ Kw-\( x )- 
As a result, we cannot use the adiabatic linear metric perturbatiorLas the metric 
of the far-zone expansion for the MiSaTaQuWa self-force in general.© However, as 
long as \(d/dX)'y(X)\ is small, the result obtained with the MiSaTaQuWa self- force is 
approximately correct. We conjecture that the condition \(d/dX)'y(X)\ < 0(1) is less 
severe than (2-8) for the usual metric perturbation scheme and the adiabatic linear 
metric perturbation is valid in a longer time scale, as we discuss in the following 
subsections. 

An adiabatic linear metric perturbation does not satisfy the linearized Einstein 
equation (2-5). Hence, it is not defined by the usual metric perturbation scheme. 
Here we briefly discuss a systematic method to construct a post-adiabatic nonlin- 
ear metric perturbation. In addition to the small parameter used in the far-zone 
expansion m/L, we introduce another small parameter for the adiabatic expansion 
e ~ |(d/dA)7| < O(l), and we expand the metric and the stress-energy tensor as 

+(m/L)h^ + (m/L)eh^ + (m/L)e 2 h^ + ■■■ 
+(m/LfhW + [m/Lfeh^ + (m/L) 2 e 2 h^ + ■■■ 

+ ••■, (3-2) 
= (m/L)T^'°^ v + (m/L^T^ 1 ^ + [m/L)e 2 T^ 2 ^ u + ■■■ 
+(m/L) 2 T^ 2 '°^ u + {m/LfeT 1 - 2 ^ + (m/ L) 2 e 2 T (2 ' 2 ^ v + ■■■ 

+ ••■, (3-3) 

In the limiting case e — > 0, this metric perturbation scheme becomes the usual metric 
perturbation scheme, and therefore some results for the usual metric perturbation 
scheme may apply to h<$ 0) and T^ n ^ u . We define to be hffi l) and T^^ v 

to be the stress-energy tensor of a monopole particle. In contrast to the usual 
metric perturbation scheme, T^' ^ u does not necessarily satisfy the conservation 
law in the background metric by itself, because the violation of the conservation law 
for r^- )^ can be compensated for by a higher order expansion with respect to e. 
Hence, the orbit of the particle described by T^'°^ v is not necessarily a geodesic of 
the background, and the expansion of the orbit in the form (2-7) is not necessary in 
this scheme. 

The violation of the conservation law presents a difficulty with regard to deriving 
equations for post-adiabatic linear metric perturbations, i.e., h$u with m ^ 1. 
Schematically, an adiabatic linear metric perturbation satisfies the equation 

G^ v [{m/L)h ad M] = (m/L)T^>°^ v + eA {1 '°^ u [(m/L)h ad ^} , (3-4) 
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where y^ 1,1 )^ appears because T^ ,0 ' liV does not satisfy the conservation law, while 
the LHS of (3-4) is algebraically divergence free. One may expect that the relation 
G^[(rn/L)eh^] = {m/L)eT^ l ^ v - {m/L)eA^'°^ u [(m/L)h ad ^} holds for the 
first post-adiabatic linear metric perturbation h^ 1,l \ where (m/L)eT^ 1 ' 1 ^ fJjU can be 
derived appropriately with a matched asymptotic expansion. However, this equation 
might be inconsistent, because the RHS is not necessarily divergence free, while the 
LHS is algebraically divergence free. We suppose that, as in the case of (3-4), there 
is a rank- two symmetric tensor that is 0((m/ L)e 2 ) and cancels the divergence of 
the RHS. As a result, we obtain the following equations for the post-adiabatic linear 
metric perturbations: 

G^ u [{m/L)e n h^] = (m/L)e n T^ n ^ u 

-eA il ' n - 1)flv [(m/L)e n - 1 h {1 ' n - 1) ] + eA^ n ^ u [(m/ L)e n h^} . (3-5) 

We believe that equations for nonlinear metric perturbations would be derived in a 
similar way. 

3.2. Adiabatic extension of the radiation reaction formula 

In order to see that the new metric perturbation scheme is valid over a longer 
time scale than the usual metric perturbation scheme, we extend the argument of 
the radiation reaction formula. In the radiation reaction formula, a component of 
a self-force averaged over an infinite time is derived in the usual metric perturbation 
scheme from the general form of the self-force. It is shown that the loss of orbital 
energy and angular momentum due to the self-force averaged over an infinite amount 
of time is equal to that by the momentum flux caused by gravitational waves. 
Since the self-force is entirely gauge dependent in the usual metric perturbation 
scheme, we can choose a gauge condition such that the rest of the self-force is small. 
Then we can say that the self-force of the radiation reaction formula is defined by the 
averaged energetic gauge condition (2-16) over a time scale for which the usual metric 
perturbation scheme is validly We consider application of the adiabatic extension of 
the linear metric perturbation (3-1) in this gauge condition. Because of the averaged 
energetic gauge condition (2-16), the radiation energy (2-12) does not decrease by 
radiation reaction, and therefore we expect that the adiabatic metric perturbation 
scheme is valid over a time scale longer than that for the usual metric perturbation 
scheme. 

Using the result given in Appendix |A|, the stress-energy tensor of a geodesic 
7 = {£ a , X b , C c } can be defined as 

T$(x) = J d\s al3 {X)5{t - J)6(r - z r )8{9 - /)<5(c/> - z*) , (3-6) 
8*0 (\) = s a ^™) expKxr + in eX e] , (3-7) 

n r ,ng 

where the expansion coefficients s a P( n r, n e) are functions of £ a . Because a Kerr black 
hole is a stationary and axisymmetric solution, we can assume that a Green function 

*' A more complete argument for this gauge condition is given in Sec.0 (see (4-11)). 
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for a linear metric perturbation in the Boyer-Lindquist coordinates can be defined 

as 

Gn VOl >{3'(x,x') = ^g^)p,{r : e]r\e')eKY>[-iuj{t-t') + im((j)-(j)')]. (3-8) 

u>,m 

The radiation reaction formula uses a linear metric perturbation derived with the 
Green method Ell It is schematically written as 

M*5 7)= E fc (a, ' m '"- ne) (7) 

xh (u,m,nr,ne) ( £ a. ^ ^ exp [_^ t + im ^] ; (3.9) 

( 7 ) = f dX l^, m ,n T ,n g )^ exp [ icJKi _ imK ^ _ irlrXr _ ingXg ] 



2^ jjLu, m ,nr,ne)/ca\ c i(uiC t -mC^) c -i2TT(nrS7rX r +n g S7 g X e ) 

x<y(w-27r^ n , Br)n9 ) J (3-10) 
^m,n r ,n e = + n r J? r + n j?6) , (3-11) 

fl<k = , \ , Q T = . ~ r . , i?e = , ~ e , . (3- 12) 



Here, i?^, f2 r and 12^ are principal frequencies of the geodesic 7, and they are observ- 
able in the form of a spectral information concerning gravitational waves at infinity. 

Although (3-9) looks general, the derivation of the linear metric perturbation 
constrains the gauge condition, because one may still add a gauge mode which is not 
a function of t — t' or <p — <j)' to the Green function. We call this a physically reasonable 
class of gauge conditions, because one can easily read the principal frequencies from 
the resulting metric perturbation (3-9), and because it satisfies the averaged energetic 
gauge condition (2-16). 

For the adiabatic approximation of the linear metric perturbation (3-1), we re- 
gard 7 as dynamical variables. With the foliation function f{x), we have 

< (1) (*)= E ^~' ne) (7(/(*))) 

LO,m,n r ,ng 

xhfy m ' nr ' ne \£ a (f(x));r,0)exp[-iut + im<t>} . (3-13) 

The adiabatic linear metric perturbation satisfies the equation (3-4). Using the 
explicit form (3T3) of the adiabatic linear metric perturbation, it is seen that 
^(1,0)^1/ ^arf(l)j j g ii near or quadratic in (d/df)j(f) or linear in (d 2 /df 2 ) r y(f). We 
find that the adiabatic linear metric satisfies the Einstein equation to an accuracy 
of 0(m/L) when \(d/d\)j\ and \(d 2 /d\ 2 )^/\ are small. In the next subsection, we 
investigate the evolution of 7 due to the adiabatic linear metric perturbation and 
show that the assumption of the adiabaticity is consistent. 

*' Because the original radiation reaction formula uses a symmetry property of the linear metric 
perturbation itself rather than the metric perturbation derived with the Green function (3-8), we 
present a derivation of the self-force in Appendix 0. 
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3.3. Adiabatic evolution of an orbit 

Here, we study the orbital evolution caused by the adiabatic metric perturbation 
(3-13). Under the assumption that \(d/d\)^\ and \{d 2 /dX 2 )^\ are small, (3-13) ap- 
proximates the linear metric perturbation, and we can still use the general expression 
of the self- force given in Ref.|^) by simply changing 7 to 7(A). This yields 

±£ a = £ ainr ' n<>) (£ a W)eMinrXr + ^n e xe} , (3-14) 

" A n r ,n e 

where Xb = 2vri? b (£ a (A))(A + A b (A)). For the orbital evolution, RefJ|) uses the 
usual metric perturbation scheme and treats the orbital evolution by expanding the 
orbit as (2-7), which is constrained by (2-8). In the framework of the adiabatic 
metric perturbation, we cannot use (2-7). Therefore we must treat the orbit in 
a nonperturbative manner. With the idea of the adiabatic evolution of the orbit, 
we consider the evolution of the orbital "constants" instead of the evolution of the 
coordinates. 

Although Refj2|) uses the perturbative expansion of the orbit (2-7), the extrap- 
olation of its result suggests 

£ a -£g ~ 0((m/L)(t/L)), 

X b -X b ~ 0{(m/L){t/L) 2 ), C c -Cl ~ 0{(m/L){t/L) 2 ) , (3-15) 

where {£$, Xq,Cq} are the initial values. Thus, the evolution of the primary "con- 
stants" £ a becomes nonperturbative after the radiation reaction time, T ra( i ~ 0(L(m/ , L)~ 1 ), 
whereas the evolution of the secondary "constants" (X b ,C c ) becomes nonperturba- 
tive after the dephasing time, Td ep ~ 0(L(m/ ' L)~ l l 2 )(< T ra d). In Appendix |C|, we 
derive evolution equations for the secondary "constants", assuming that the evolu- 
tion of the primary "constants" is perturbative, and we find that the approximation 
(3-15) obtained with the usual metric perturbation scheme is qualitatively correct. 

Equation (3-15) shows that the assumption of the adiabaticity is constrainted by 
the conditions \(d/d\)(\ b , C c )\ < O(l) and \(d 2 /d\ 2 )(X b , C c )\ < O(l), and the adi- 
abaticity holds on the radiation reaction time scale, where the primary "constants" 
still evolve perturbatively. Hence the adiabatic linear metric perturbation satisfies 
the Einstein equation to an accuracy of 0(m/L) over the radiation reaction time 
scale t < T ra( j. 

In Appendix O we integrate the evolution equations in the case that the orbit is 
either eccentric in the equatorial plane (9 = vr/2) or circular in an inclined precessing 
plane (r = const). In either case, it is shown that the circular (equatorial) orbit 
remains circular (equatorial) under the influence of the self- force. In Appendix 
, we find that the self-force in the adiabatic linear metric perturbation scheme can 
be used to predict the orbit over the radiation reaction time scale. The dominant 
part of the orbit is described by only the zero mode of the self-force, £ a (°'°) ; which 
is the gauge invariant component of the self-force in the usual metric perturbation 
scheme, a) This is a reasonable result of the adiabatic metric perturbation scheme. 
Because the adiabatic linear metric perturbation (3-13) is a valid metric perturbation 



Self-Force in the Radiation Reaction Formula 



15 



on this time scale, the background metric accurately approximates the full geometry, 
and the orbit using the background coordinates as an approximate reference must 
be gauge invariant to leading order. We also find that this is observable. Because 
the adiabatic metric perturbation scheme describes the metric perturbation together 
with the orbital evolution caused by the self-force, the modulation of gravitational 
waves induced by the radiation reaction can be read off of (3-13). In (3-13), we have 

k (uJ,m,n r ,ng)^^ ^ e -i2Tr{n r Q r \ r (f)+n e Sl e X s (/)) ^(wC* (f)-mC* (/)) ^ /g.-^ 

and the nonperturbative evolution of (X b ,C c ) can be observed as a modulation of 
the wave phase. 

The remaining component of the orbit involves nonzero modes of the self-force 
(3-14), £ a K> n e), with (n r ,ne) ^ (0,0). They impart an 0((m/L)) contribution to 
the orbit on the radiation reaction time scale. This is the same order as the metric 
perturbation, and therefore these modes are interpreted as a gauge. In the next sec- 
tion, we show that the nonzero mode of the self-force £ a ( nr > n e) ; with (n r ,ng) ^ (0, 0) 
is gauge dependent in the adiabatic metric perturbation scheme. We also find that 
the adiabatic nonlinear metric perturbation imparts an 0((m/L) 2 (t/L) 2 ) contribu- 
tion to the orbit. Well within the radiation reaction time scale, this is smaller than 
O(l), and it is consistent with the fact that the adiabatic linear metric perturbation 
satisfies the Einstein equation to an accuracy of 0(m/L) on the radiation reaction 
time scale. 

Although the above may seem quite reasonable, there is an additional concern 
when we consider a general orbit. By looking at (3-14), one can easily understand the 
reason why the gauge invariant component of the self- force £ a (°'°) imparts a dominant 
contribution to the orbit. In the adiabatic metric perturbation scheme, £ a (°'°) slowly 
changes over the radiation reaction time scale and it coherently contributes to the 
orbital evolution. Contrastingly, the nonzero modes of the self- force, (n r ,ng) ^ (0, 0) 
give a small contribution to the orbital evolution, because the phase n r \ r + ngxe 
changes on the dynamical time scale in general, and the contribution cancels over 
a longer time scale. This is true for a circular or equatorial orbit. However, there 
could be a case in which the phase n r x r + nexe of a nonzero mode becomes almost 
stationary. 

This problem appears as a discrepancy between the gauge invariant quantity 
in the radiation reaction formula and in the balance formula in the framework of 
the usual metric perturbation scheme. In Refj2[), the gauge invariant component 
of the self-force is derived by taking an infinite time average over the self-force. 
The infinite time averaged self-force is derived in the usual metric perturbation as 
-^£ a ^j = £ a ( Ur > n0 \ where the sum is taken over the values of n r and uq such that 

n T Q T + ngfig = 0. On the other hand, Ref.[l(]) shows that the zero mode, £ a (°> \ 
is equal to the gauge invariant momentum flux of gravitational waves. In order to 
resolve this disagreement, it is assumed in Refj2|) that Q r / Qq is irrational. In this 
case, the infinite time averaged component of the self-force becomes equal to the 
gauge invariant momentum flux of gravitational waves. However, this assumption 
is not strictly valid. If there is a nonzero mode for which n r Q r + ngilg = 0, the 
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phase n r Xr + n eXe becomes almost stationary, and it may contribute to the orbital 
evolution. In this case, we may need to develop another formula, since we can only 
calculate £ ,a ( 0,0 ) using the radiation reaction formula. 

By extending the analysis of Appendic |D|, we conjecture that such a problem 
might not arise in the adiabatic metric perturbation scheme and that all we need for 
an orbital prediction is £ a (°'°). The primary "constants" for a general orbit during 
the radiation reaction time, t < 0(L(m/ X) -1 ), evolve according to 

£ a = £$ + [ X d\ £ a(nr ' ne) (£o) expKxr + inexe] + 0(fi 2 t) . (3-17) 

*'^0 n r ,n g 

During the dephasing time interval < t < 0(L(m/ X) -1 / 2 ), we have Xb ~ + 
const. As a result, the phase of any nonzero mode with | n r J7 r +ne/?e | < 0(1/ 'L(m/ X) 1 / 2 ) 
becomes stationary and gives a coherent 0((m/ X) 1 / 2 ) contribution to the primary 
"constants" given in (3-17). Beyond the dephasing time, Xb begins to evolve in a 
nonperturbative manner, and, even if \n r fi r + riQf2g\ is small, n r x r + n eXe would 
not be stationary unless (n r ,n$) = (0,0). During the time interval satisfying 
0(L{m/L)- 1 / 2 ) < t < 0(L(m/L)- 2 / 3 ), a mode for which \{n r Q^ a +n e Q e ,a)£ a{afi) \ < 
0(1/ L(m/ X) 2 ) would have a stationary phase. Extending this approximation, dur- 
ing the time interval 0(L(m/L)~ n ^ n+1 ^) < t < 0(L(m/L)-( n+1 )/( n+2 )), a mode 
for which \(n r n r>a + n e ^ ia )(^ a (°'°)) n | < 0(1/L{m/L) n+1 ) would have a stationary 
phase, and it would contribute an 0((m/X) 1 /("'+ 2 )) < O(l) term to (3-17) during this 
time interval. Thus, only the zero mode £ a (°'°) could impart a coherent contribution 
of 0(1) over the radiation reaction time scale. 

Before further investigating the gauge condition, we discuss the validity of the 
adiabatic metric perturbation for an application to the LISA project. The LISA 
project is planning to detect gravitational waves for several years, and we need 
theoretical templates of gravitational waveforms valid over such time intervals. It is 
estimated that the mass ratio of a possible target binary would be ~ 10 -5 . Suppose 
m/L = 10~ 5 with a dynamical time scale L = 100s. Then the radiation reaction 
time T ra( i would be on the order of several months. Hence, a theoretical prediction 
obtained using the adiabatic linear metric perturbation may not be reliable enough 
for the LISA project. 

In the next section, we argue that fixing the gauge condition improves the va- 
lidity of an adiabatic linear metric perturbation, and with this improvement one can 
calculate gravitational waveforms with the adiabatic metric perturbation in a more 
reliable manner. 

§4. Radiation reaction gauge 

In this section, we consider a further gauge transformation of an adiabatic met- 
ric perturbation. It is found in Ref.|2|) that in the usual metric perturbation scheme, 
the gauge transformation of the self-force acting on the primary "constants" can be 
written in terms of a total derivative along the orbit. For this reason, the gauge 
transformation of the primary "constants" can be expressed locally. Because the 
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argument used here depends on only the properties of the orbital equation, the 
same result can be used for an adiabatic metric perturbation. However, for techni- 
cal simplicity, here we first consider the gauge transformation in the usual metric 
perturbation scheme. Then we extend the result to the case of an adiabatic metric 
perturbation. 

Suppose we transform the coordinates as — » x M + Then, to leading order 
in the gauge vector we have the primary "constants" transform as 

8S E / L = -rfJ L v^% + vZ L v^ a , 5C C = -rfc^t" , (4-1) 

TP /T 

where v a is the 4- vector of the orbit, rj a are the Killing vectors used to define the 
orbital energy and the z component of the angular momentum, and rj^ is the Killing 
tensor used to define Carter constant. In the usual metric perturbation scheme, the 
primary "constants" evolve according to the relation 

£ a = £g + £ a{m X 

+ E , 7 1 ^£ a(nr > ne) exp[i(n r n r + in e n e )\] , (4-2) 

(nr.n^co.o) i{n r n r + n e ne) 

where we have excluded the exceptional case that Q r /Qg is rational, as in Ref.^). 
We now consider the elimination of the last oscillating part of (4-2) through a gauge 
transformation. Although £ Q is a vector field defined in the entire spacetime, the 
following analysis considers the condition for £ a along the orbit. Therefore, we 
consider £ a to be a function of the orbital parameter as £ a = £ Q (z(A)). 
As is shown in Appendix ^, the gauge equation (4-1) becomes 

d 
'dX 



b=E,L,C 



5£ a = J2 (A a b —e + B a b e), (4-3) 



where A is a regular 3 x 3-symmetric matrix and B is a 3 x 3 anti-symmetric matrix. 
With these properties, one can formally integrate the gauge equation as 



£ = CT 1 / d\CA~ 1 5£ , (4-4) 
where C is a matrix schematically defined by 

C = exp( J dAA _1 B) . (4-5) 

We recall that, because the components of the matrix A _1 B are functions of r(A) 
and 0(A), their components can be Fourier-decomposed as exp(in r Xr + ingxe) (see 
Appendix |A]). Applying an appropriate mathematical transformation, the matrix C 
becomes 

C = exp(iciU! + ic 2 u 2 ) , (4-6) 

where ui and u 2 are two independent elements of the so(3)-algebra, and two coeffi- 
cients ci and c 2 are functions of A that satisfy 

ci/2 = ci/ 2 A + ^2 ci/2' n<,) exp(m r Xr + inexe) ■ (4-7) 
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For this reason, the matrix C can be decomposed into discrete Fourier modes, and 
we schematically write 

C = ^C w exp(iu;A). (4-8) 

w 

In order to eliminate the oscillating term of (4-2), we define 5£ a for (4-4) as 
5£ a = £a(nr,ne) expKxr + inexe] - 5£g , (4-9) 

n r ,ng 

where 5£§ is a constant. If we set 8£§ = 0, we would have £ a ~ 0((m/L)X), due 
to the integration over A in (4-4), and as a result, the gauge transformation of the 
metric perturbation would grow as 0((m/L)(t/L)) around the orbit. Such a gauge 
transformation would cause the metric perturbation to become invalid in a short 
time in the adiabatic metric perturbation scheme. In order to avoid this, we choose 
5£q 7^ 0, so that the gauge transformation (4-4) becomes oscillatory. Then we have 

£ a = {£$ + 5£g) + £ a( -°^X . (4-10) 

This shows that nonzero modes, £ a ( n r> n e) with (n r ,ng) 7^ (0, 0), are completely gauge 
dependent. We recall that the initial value £q has a 0(m/L) ambiguity when we 
include the metric perturbation. Therefore, 5£q can be renormalized into the initial 
value, £q . 

We now consider the self-force in the adiabatic metric perturbation scheme. The 
self-force under this gauge condition becomes 

d £a = ^(0,0) _ (4 . n) 

dX 

We refer to this gauge condition as the "radiation reaction gauge", as the self-force 
is written solely in terms of the gauge invariant dissipative term. IUs interesting to 
note that a previous calculation employing the balance formulaLp'EP uses the same 
equation. However, it was considered to be an approximation and a validity of this 
approximation has been discussed. 0) The derivation here uses the gauge freedom 
of the metric perturbation. Hence, it is an exact result in the framework of the 
adiabatic metric perturbation. 

The application of this result to a motion of a spinning particle is interesting. A 
spinning particle does not move along a geodesic due to the spin-curvature coupling, 
t3 and one may regard this effect as a force: 



D 

dX 



v a = F? pin . (4-12) 



However, this is a conservative effect, and the conservative force can be eliminated by 
choosing the radiation reaction gauge, as long as the force is on the order of the metric 
perturbation, F^ pin ~ 0{m/L). This is the case for a spinning compact object, and 
the result suggests that we may not observe the spin effect by gravitational wave 
detection. Although the force resulting from the spin-curvature coupling can be 



Self-Force in the Radiation Reaction Formula 



19 



eliminated, this does not mean that there is no such effect. As (4-10) shows, it is 
simply absorbed into the initial value. As a result, it alters the condition for the last 
stable orbit of the binary. 

It is interesting that the radiation reaction gauge changes the validity argument 
for an adiabatic linear metric perturbation. Because {d/dX)£ a does not have nonzero 
modes, we have {d 2 / dX 2 )£ a ~ 0((m/ L) 2 ) under this gauge condition. We also have 
((d 2 /d\ 2 )X b , {d 2 / d\ 2 )C c ) ~ 0(m/L) and the validity of the adiabatic linear metric 
perturbation depends on only ((d/dX)X b , (d/d\)C c ). 

Here we use the post- Newtonian approximation for (d/dX)£ a = £ a (°'°) ^ be- 
cause £ a (°> ) is a throughly studied gauge invariant quantity. Although the post- 
Newtonian expansion is not a good approximation for a strongly gravitating system, 
such as an extreme mass ratio binary system, an extrapolation is believed to give 
a good approximation of the result. In the perturbation limit, we have £ a (°'°) ~ 
0((m/ L)v 5 ), where v is a typical velocity of the system, and v 2 is estimated to be 
less than 0.1 during the inspiral stage of the binary. We have ((d/dX)X b , (d/dX)C c ) ~ 
0((m/L)v 5 (t/L)), and therefore the assumption of adiabaticity holds for the post- 
Newtonian radiation reaction time scale, t < TpN ra( i ~ 0(Lv~ 5 (m/ 'L) _1 ), over 
which the adiabatic linear metric perturbation could be a solution of the Einstein 
equation to an accuracy of 0(m/L). 

As in the previous section, we suppose m/L = 10~ 5 with a dynamical time scale 
L = 100s for a possible target of the LISA project, and we have Tpjv.md on the order 
of several years. Hence, we believe that the adiabatic linear metric perturbation 
in the radiation reaction gauge should give us a reliable prediction of gravitational 
waveforms for the LISA project. 

§5. Summary 

In this paper, we have considered the so-called self-force problem in the calcula- 
tion of gravitational waves from an extreme mass ratio binary. A strategy previously 
employed for this problem was that, because we already have a method to calculate 
gravitational waves for a given stress-energy tensor, we calculate a self-force in order 
to derive the orbit of a particle. A general form of a self-force called the MiSa- 
TaQuWa self-force, was derived in the usual metric perturbation scheme. a> We find 
that, in the usual metric perturbation scheme, the self-force is completely gauge 
dependent, and it can even be the case that the self-force vanishes. 

In Sec|2], we argue that this is because the usual metric perturbation scheme 
is defined to allow an arbitrary gauge condition, and as a result, an inappropriate 
choice of the gauge condition would eventually lead to a nonperturbative evolution 
of the metric from an energetic point of view. Such a consequence cannot be avoided 
in the usual metric perturbation scheme, because the scheme requires that the linear 
metric perturbation must be induced by a background geodesic, and therefore the 
metric perturbation cannot be tracked with the orbital evolution caused by the self- 
force. Although the nonperturbative evolution of the metric could be calculated using 
numerical relativity, it would require a very large computational cost, and for this 
reason, we consider it advantageous to modify the metric perturbation scheme with 
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a constraint of a gauge condition for the self-force. In Sec||, we propose an adiabatic 
metric perturbation scheme. In this new metric perturbation scheme, the so-called 
adiabatic approximation is applied to both the metric perturbation and the orbit. As 
a result, the metric perturbation is derived together with the orbital evolution caused 
by the gravitational self-force. We find that the new metric perturbation scheme is 
valid over the radiation reaction time scale, whereas the usual metric perturbation 
scheme is valid over only the dephasing time scale, which is shorter than the radiation 
reaction time scale by the square root of the mass ratio (~ 10~~ 2 ' 5 ). 

The adiabatic approximation has been successfully applied to linear metric per- 
turbations, but it is not a simple problem to apply it to post-adiabatic metric per- 
turbations because the linearized Einstein operator is algebraically divergence free, 
as we discussed briefly in Sec||. Although it is valid for sufficiently small (d/dX)j 
and (d 2 /dX 2 )>y, this may pose a question regarding the convergence of the adiabatic 
extension. 

For the self- force problem, there are two approaches: One employs the self- force 
calculation, which involves an explicit regularization calculation, and the other 
employs the radiation reaction formula,^ which is a nontrivial extension of the 
balance formula. & In this paper, we have defined the adiabatic metric perturbation 
scheme using the radiation reaction formula. One may ask whether an adiabatic 
extension is possible for the self-force calculation. We believe that this would be 
nontrivial, because, for a technical reason, the self-force calculation cannot be easily 
carried out with a linear metric perturbation taking the form (3-9). When the 
background is a Schwarzschild black hole, one may use the Zerilli-Regge- Wheeler 
formalism© to calculate a metric perturbation. In this case, one can transform the 
metric perturbation into the form (3-9). However, when the background is a Kerr 
black hole, a method to calculate an inhomogeneous metric perturbation has only 
been proposed^ and it has not yet been implemented numerically. It is not clear 
that the metric perturbation in this method can be transformed into the form (3-9). 

Because all we need for the orbital evolution in the end is a gauge invariant 
quantity £ a (°>°) , one may consider the possibility of calculating it by taking a long 
time average of the result of the self-force calculation, & assuming that the self- force 
calculation is implemented under an appropriate gauge condition. Although this 
may sound a plausible approach, the self-force averaged over a long time T becomes 



dX I x 



where the sum is taken over n r and rig for which \n r fl r + ngf2g\ < 1/T (see (3-14)). 
Thus, there is a case in which this approach fails to provide £ a (°>°). 

By contrast, the radiation reaction formula^ requires only a homogeneous met- 
ric perturbation and a conventional method for this calculation is known. cZP Because 
the necessary calculation is only a slight extension of the balance formula, we con- 
sider the numerical technique for this calculation to be established. A homogeneous 
metric perturbation can be derived by operating with a tensor differential operator 
on the homogeneous Teukolsky function, HZP and an analytic method for the homoge- 
neous Teukolsky function is known. Because of its fast convergence, this analytic 
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method gives a more accurate result in a more efficient way. US' 

In Sec.|], we introduced the radiation reaction gauge, in which the self-force is 
expressed exactly by the radiation reaction formula (4-11). We recall that the ener- 
getic argument given in Sec.|2] shows that the metric could be nonperturbative with 
an inappropriate choice of the gauge condition for the self-force. Because the usual 
metric perturbation scheme allows an arbitrary gauge condition, it is valid only over 
a short time scale. Contrastingly, the validity of the adiabatic metric perturbation 
depends on the gauge condition, because it solves the metric perturbation and the 
orbit at the same time. We find that the adiabatic metric perturbation under the 
radiation reaction gauge condition is valid over a longer time scale than that under 
a general gauge condition for the radiation reaction formula (3-9). We also note that 
the evaluation of the validity might not be correct for a highly eccentric orbit. 

It is also interesting that a previous study using the balance formula lP'EP uses the 
equation (4-11) by the assumption of the adiabatic approximation. In this approach, 
the adiabatic approximation holds if the orbital period is sufficiently shorter than 
the dephasing time. However, the analysis in Sec.|| shows that one can consistently 
eliminate the oscillating part of the self-force through the gauge, as long as it is 
small. Hence, the condition for (4-11) must be the same for the adiabatic metric 
perturbation scheme. 
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Appendix A 

— Geodesic — 



We consider inspiralling geodesies in a Kerr geometryH A geodesic equation is 
a set of 2nd order differential equations. Therefore, geodesies are characterized by 
6 integral constants. 3 of them are constants of motion, the orbital energy (E), the 
z-component of the angular momentum (L) and the Carter constant (C). In order 
to distinguish the rest of the constants, we call E, L and C the primary constants 
and denote them by £ a , a = {E,L,C} when it is appropriate. Using the primary 
constants, the geodesic equation becomes 



dr\ 2 



- 



dXJ 



(E{r 2 + a 2 ) - aL) 2 - (r 2 + (L — aE) 2 + C)A , (A-l) 



*' A rather sophisticated analysis using Hamilton- Jacobi formalism is given Refjl9[). Here we 
summarize another derivation given in Refj2|) because it is convenient for the analysis given in this 
paper. 
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C-a 2 (l-E 2 )cos 2 e-L 2 cot 2 9, (A-2) 
-a(aE sin 2 - L) + (r 2 + a 2 )(£(r 2 + a 2 ) - aL)/Z\ , (A-3) 

-{aE - Lcosec 2 <9) + (a/A)(E(r 2 + a 2 ) - aL) , (A-4) 
r 2 - 2Mr + a 2 . (A-5) 

Here we use A as an orbital parameter which has a one-to-one relation with the 
proper time, r, given by r = / e?A(r 2 + a 2 cos 2 #) 2 . 

For inspiraling geodesies, r/#-motion is bounded and the solutions of (AT) and 
(A-2) can be expanded in discrete Fourier series as 

z b = j2 z b ^ exp[m X 6] , Xb = ^Q b {\ + X b ) , (A-6) 

n 

where b = {r, 9}, with r = z r and 8 = z e , and Z b ^ and 1?^ are functions of £ a . Here 
we have two integral constants X b , but, since we can freely choose the zero point of 
A, A r — X e alone specifies the geodesic. (A-3) and (A-4) are easily integrated, and we 
obtain 

z c = k c + J2 Z b {n) e MinXb] , k c =<Z c >X + C c , (A- 7) 

b,n 

where c = {t, <p} with t = z t and (j> = z^, and (z c ^j and are functions of £ a . 

Here we have two integral constants C c which specify the geodesic. We call X b and 
C c the secondary constants of the geodesic. 

Because of the periodicity, one can freely add 1/17& to X b . In general, the ratio 
of tl T to fig is irrational, and one can take A r — A as small as possible using this 
freedom. The constants C c can be set zero by using the t/^-translation symmetry 
of the Kerr geometry. 

Appendix B 

Self-Force 

We argue that a metric perturbation defined by (3-9) yields a general form of 
the self-force given in Ref.|2|). In Refj2|), a general form of the self-force (3-14) is 
derived, using a symmetry property of the Kerr geometry and a symmetry property 
of a geodesic. Because the self-force is induced by a geodesic, it is a function of the 
geodesic parameter A and the geodesic constants 7 = {£ a ,X b ,C c } (see Appendix 
[A| for definitions). Employing the symmetry of the Kerr geometry, Ref.^) finds 
that the self-force can be independent of C c . Then by a symmetry property of a 
geodesic, the self-force can be expanded in a discrete Fourier series of terms of the 
form exp[n r (A — A r ) + ng(X — X e )]. In this appendix, we demonstrate that the gauge 
condition for the metric perturbation (3-9) actually respects the symmetry of the 
Kerr geometry by showing that the C c derivative of the self- force vanishes. 



d9\ 2 
dXJ 
dt 
dX 



dX 
A 
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We assume that the self-force acting on the orbital energy, angular momentum 
and Carter constant can be derived using a regularization calculation Ef as 

^-£ a (X)= Hm F>(x; 7 )-H,; 7 )], (B-l) 

where -F a [/i] is a tensor derivative operator acting on the metric perturbation h^ u . 

and z^(A;7) are a field point and an orbital point of the geodesic 7 respectively. 
h^ u (x;j) is defined by (3-9). h s ™ 9 (x;j) is called the "singular part" of the metric 
perturbation, and it is defined in Ref.[|) under the harmonic gauge condition. 

Using the temporal and rotational Killing vectors c = {t, (/>}, the symmetry 
property of the Kerr geometry is defined in terms of the Lie derivative along the 
Killing vectors as 

%)flSr = 0. (B-2) 

The property (3-8) of the full Green function for a metric perturbation is expressed 
as 

(£s {c) +£'s (c) )G^ af i(x,z)=0, (B-3) 

where £f. . acts on x and the indices \i and u, and £'§ (c) acts on z and the indices 

a and 0. The property (3-8) is consistent with the linearized Einstein operator G$ 
since it commutes with the Killing vectors, i.e. 

[£ kc) ,G$] = 0. (B-4) 

It is possible to derive h s ™ 9 (x\ 7) using the Green method with a singular Green 
function, and an explicit form of the singular Green function in the harmonic gauge 
condition is constructed using a bi-tensor expansion in Ref.[|). Because bi-tensors 
are defined in terms of background geometric quantities, the singular Green function 
in the harmonic gauge condition can also satisfy the relation 

(A( c) +^ (c ,)G^(x,z)=0. (B-5) 

As in Re£pO|), we assume that a gauge transformation acts only on the full Green 
function.fcjThen, the singular Green function has the property (B-5) in general. 

Because the singular behavior of the full Green function cancels with that of the 
singular Green function, one can define a regular Green function by 

Gjfttffa z ) = G^a/3(x, z) - G s ™^(x, z) , (B-6) 



''■-Strictly speaking, there are two types of gauge transformations for the MiSaTaQuWa self- 
force. EP One is for the metric of the far-zone expansion and is discussed in Refj2C|). The other is for 
the metric of the near-zone expansion. The metric of the near-zone expansion is defined as the sum 
of the black hole metric as a point particle and its perturbation. In RefjE]), we choose the gauge 
condition for which the particle is located at r = of the Schwarzschild coordinates when the L — 1 
mode of the black hole perturbation vanishes. With this condition, the singular part is uniquely 
determined. 
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and the self-force is simply derived as 

±£ a (\) = F a [h re °](z(\y n ) , (B-7) 
h™ 9 (x;j) = J dz^ afi (x,z)T^{z-i) , (B-8) 

where T a/3 is the usual stress-energy tensor of a point particle moving along the 
geodesic 7. 

Using (B-3) and (B-5), and integrating by parts, we have 

£ ?(c) 7) = / dz'G^ a ^z)C! kc) T^{ z -i) . (B-9) 
Prom the result of Appendix we find 



C' kc) T^(z; 1 ) = ---T^(z; 1 ). (B-10) 



Since the C c derivative of the self-force (B-7) acts on z(X) and 7 as 



x=z(A) 



= 0, (B-ll) 
we find that the self-force is independent of C c . 

Appendix C 

Evolution Equations 



The evolution equations for the primary "constants" are given by the self-force 
(3-14). In this appendix, we derive the evolution equations for the secondary "con- 
stants" (A b ,C c ). For convenience, we define the orbit as 



z b = j2 exp[inx 6 ] , z c = R c + ]T zf n) exp[mx&] , (C-l) 



n b,n 

where Z b ^ and Z 6 c(n) are the functions of £ a given in (A-6) and (A-7), and consider 
the evolution equations for (xb, & c ). We can derive the evolution for (A 6 , C c ) by using 
Xb = 2Trf2 b (X + X b ) and k c =< Z c > A + C c (see (A-6) and (A-7)). 
Using (xb,K c ), the self-force (3-14) is written as 

±£« = F a = £ a{nr ' ne) (£ a )eMin r Xr + ineXe] - (C-2) 

dA n r ,n g 

It is important to note that (A-l)-(A-4) hold for non-geodesic orbits, because they 
are derived from the definitions of (E, L, C), whether (E, L, C) are constant or not. 
Together with (C-2), (A-l) and (A-2) give the evolution equations for Xbi and we 
obtain the evolution equations for k c from (A-3) and (A-4). 
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We rewrite (A-l) and (A-2) as 

{iff < c ' 3 > 

For convenience, we define 

5 x z b = in2nf2 b Z b ^ exp[in X b] , 6 a z b = £ Z$W exp[inx b ] , (C-4) 

n n 

where j0 represents the derivative with respect to £ a . Then, we have 

dz h 5 x z b fdx b 



(f?) +£(^ a , (C-5) 



dX 27rl?b \ d\ 
and, from (C-3), we obtain the evolution equations for Xb as 

where we have used (S\z b ) 2 = Z b (£,z b ). 

We next consider the t/4>- motion. We rewrite (A-3) and (A-4) as 

% = zZm,z h )- (C-7) 

For convenience, we define 

5 6 z c = J] m27r4^ c(n) exp[mx 6 ] , <U C = £ exp[inx 6 ] ■ (C-8) 

n b ,n 

Then, we have 

and, from (C-7), we obtain the evolution equations for K/£ clS 

*= W + ? (^. + ? *0!)#., (,..) 

where we have used + Sft = J2b -^b(£; z b ). 

In summary, we obtain the evolution equations for (X b , C c ) as 

It is notable that, although we have the freedom to choose the zero point for A and 
to add 1/ ,17& to X b , the last terms of (CT1) and (CT2) eventually grow linearly in A. 
Therefore, ((d/dX)X b , (d/dX)C c ) are estimated to be 0((m/L)(t/ L)), and similarly 
((d 2 /dX 2 )X b , (d 2 /dX 2 )C c ) behaves as 0((m/L)(t/L)), 
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Appendix D 

Orbital Evolution: Circular or Equatorial Case 

We consider integration of the orbital equations (C-2), (C-6) and (CTO) in order 
to understand the properties of orbits. Here we consider only either an eccentric 
orbit in the equatorial plane (0 = vr/2) or a circular orbit in an inclined precessing 
plane (r = const). In this case, it is sufficient to consider just one of Xb (or At). 
We also use the fact that the primary "constants" can be treated in a perturbative 
manner over the radiation reaction time scale. 

With an appropriate scaling, the evolution equations (C-2), (C-6) and (C-10) 
can be reduced to 

^ = A + 0(a 2 ), A = Y i a n e in \ (D-l) 
^ = 1 + x + B + 0( Ji 2 ), B = J^b n e m y, (D-2) 
= l + x + C + 0(i?), C = J2cne my , (D-3) 

at n 

where (x, y, z; t) corresponds to (£ a , x b , H c ; A). Here, n = m/L < 1 is used as an index 
of the approximation. Thus, (a n ,b n ,c n ) ~ O(ji). The 0(/j, 2 ) terms of (D-l)~(D-3) 
come from the nonlinear metric perturbation. 

For simplicity, we assume the initial values x = y = z = 0&tt = 0. From (D-l) 
and (D-2), we have 

dx , . ^, 2 



and therefore we have 



= A + 0(^), (D-4) 
dy 



x = AW + 0(n 2 y) , (D-5) 



where 

AW =a y + J2 a n ]einy 



l 1] = -E al n ] . = ^ - (n + 0) (D-6) 



in 

From (D-3) and (D-2), we have 

dz 

_ = 1 _ B + C + 0(u 2 y), (D-7) 

and therefore we have 

z = y- BW + CW + 0(/x 2 y 2 ) , (D-8) 
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where 



my 



b [1] 
°0 



Jl] _ 



n 

-E4 11 , 4 1] = ^c n . (n^O) 

n^O 



(D-9) 



(D-10) 



(D-5) and (D-8) show that the evolution of x and z with respect to y is dominatly 
determined by ao,6o and Co when y > 1 and y > fx, respectively, and the effect of 
a n , b n and c n with n ^ is only 0(/x). We also note that the effect of the self-force 
by a nonlinear metric perturbation becomes 0(1) when y ~ and therefore the 
approximation loses validity. 

(D-2) becomes 



(l-AW-S)| = l + 0(^), 
and we can integrate to obtain 

y_ A w_ B m = t + o( fl h 2 ), 

where we define A® and as 



(D-ll) 
(D-12) 



A [2] 



aoy 



/2 + 4 1] y + E«L 2] ^, 



laW(n^O), 4 21 = -E 



J 2 ! 



(D-13) 



n^O 



We can rewrite (D-12) as 

i - 4 1] " b o 
y = 

a 



2a 



(1 - of ] - 6 ) 



-(t + D) 



+ 0( M V), (D-14) 



where we define D = J2 n ( a ^ + &n')e m!/ . 

(D-14) shows that the evolution of y with respect to t is dominantly determined 
by ao,a$ and &o when t > fi, and the contribution from a n ,b n and c n with n / 
is smaller by a factor of O(pft) relative to the leading order. The approximation 
becomes invalid when t ~ due to the non-linear self-force. Since y ~ 
for £ ~ /Lt -1 , the linear self- force can be used to obtain a prediction of the orbital 
evolution for the time satisfying t < /i -1 . One can further extend (D-14) as 



-(* + £>) + 



a 



(t + J D)" + --- + 0(^ 2 i 2 )(D-15) 



2(1 - af 1 - 6 ) 3 

2.3 



t + K iJ + 6 )i + aor/2 + ait A jl + ■■■ + Q(ji, \i l t l ) . 



(D-16) 
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Since aj^ depends on the choice of initial values, this shows that only ao is essen- 
tially important in (D-l). This also suggests that (3T5) is a qualitatively good 
approximation over the radiation reaction time scale, because ant 2 /2 > a^i 3 /2. 

Appendix E 

Gauge Transformation 

To evaluate (4-1), we choose the gauge transformation vector as 

=$ E ri Ea + Z L ri La +£ c ri Ca f,v f> . (E-l) 

Then, the transformation of the primary constants defined in the perturbed space- 
time becomes 

5^= ]T (A\j-t? + B a b A, (E-2) 

b=E,L,C V dX J 



where the 3x3 matrices A a \ } and B a b are defined as 

/ ^E^Ea ^E^La ^E^C a „S 



A a h 



B a h 



dX 



-VaV La ~VaV La -rfcrf a pv? I ^ , (E-3) 

y -r] E r) Ca pvP -r] L ri Ca pvP -ri^rf a ^ ' ' 

/ v a [rj E ,ri L ] a rf^^rfP 

-v a [rj E ,V L ] a r j ^ /3 v a §r l L ^ j . (E-4) 

, _„C7 aD_ Ef3 _ n C v aR„L(3 n 



For a general orbit, one can assume that the symmetric matrix A has an inverse 
matrix, since the appropriate linear combinations ofn Ea ,r] La and rf a pv^ give three 
independent vectors. When the orbit is equatorial (or quasi-circular), the symmetric 
matrix could be singular. However, since it is sufficient to consider {E, L} (or {E, C}) 
in this case, the relevant 2x2 sub-matrix of A still have an inverse. 

References 

[1] S. Detweiler, E. Messaritaki and B. F. Whiting, Phys. Rev. D 67 (2003), 104016. 
L. Barack and A. Ori, Phys. Rev. Lett. 90 (2003), 111101. 
L. Barack and A. Ori, Phys. Rev. D 67 (2003), 024029. 
Y. Mino, H. Nakano and M. Sasaki, Prog. Theor. Phys. 108 (2003), 1039. 
L. Barack and C. O. Lousto, Phys. Rev. D 66 (2002), 061502. 
L. Barack and A. Ori, Phys. Rev. D 66 (2002), 084022. 

L. Barack, Y. Mino, H. Nakano, A. Ori and M. Sasaki, Phys. Rev. Lett. 88 (2002), 091101. 
[2] Y. Mino, Phys. Rev. D 67 (2003), 084027. 
[3] S. A. Teukolsky, Astrophys. J.185 (1973), 635. 

W. H. Press and S. A. Teukolsky, Astrophys. J. 185 (1973), 649. 

S. A. Teukolsky and W. H. Press, Astrophys. J. 193 (1974), 443. 
[4] The missio n description and its scien ce objectives of LISA project can be seen on the LISA 



Web page: http://lisa.jpl.nasa.gov/ 
[51 E. Poiss on. Living Rev. Relativity 7, (2004), 6. [Online article]: jittp: / / www. livingre vie ws.org/ lrr- 



2004-6/ 



Y. Mino, M. Sasaki and T. Tanaka, Prog. Theor. Phys. Suppl. 128 (1997), 373. 
Y. Mino, M. Sasaki and T. Tanaka, Phys. Rev. D 55 (1997), 3457. 
T. C. Quinn, and R. M. Wald, Phys. Rev. D 56 (1997), 3381. 



Self-Force in the Radiation Reaction Formula 



29 



[6] R. A. Isaacson, Phys. Rev. 166 (1968), 1263. 

R. A. Isaacson, Phys. Rev. 166 (1968), 1272. 
[7] For recent calculations of high precision, K. Glampedakis and D. Kennefick, Phys. Rev. D 

66 (2002), 044002. 

S. A. Hughes, Phys. Rev. D 64 (2001), 064004. 

S. A. Hughes, Phys. Rev. D 61 (2000), 084004; Erratum-ibid. D 63 (2001), 049902. 
[8] Analytic methods are summarized in; Y. Mino, M. Sasaki, M. Shibata, H. Tagoshi and 
T. Tanaka, Prog. Theore. Phys. Suppl. 128 (1997), 1. 

M. Sasaki and H. Tagoshi. Living R ev. Relativity 6, (2003), 6. [Online article]: 



http://www.livingreviews.org/lrr-2003-6/ 



[9] W. G. Dixon, in Isolated Gravitational Systems in General Relativity ed. J. Ehlers (North- 
Holland, Amsterdam, 1979), 156. 

[10] D. V. Gal'tsov, J. Phys. A 15 (1982), 3737. 

[11] A. Ori and E. Rosenthal, a talk in Capra Meeting at UTB. 

[12] D. Kennefick and A. Ori, Phys. Rev. D 53 (1996), 4319, 

Y. Mino, M. Sasaki, M. Shibata, H. Tagoshi and T. Tanaka, Prog. Theore. Phys. Suppl. 
No. 128 (1997) Chapter 1 Sec. 10. 

[13] A. Papapetrou, Proc. R. Soc. London, A 209 (1951), 243. 

[14] L. Blanc het. Living Rev. Relativity 5, (2002), 3. [Online article]: fittp: / /www.livingreviews.org/lrr^ 



2002-3/ 



and the reference therein. 
[15] F. J. Zerilli, Phys. Rev. D 2 (1970), 2141. 
[16] A. Ori, Phys. Rev. D 67 (2003), 124010. 
[17] P. L. Chrzanowski, Phys. Rev. D 11 (1975), 2042. 
[18] R. Fujita and H. Tagoshi, Prog. Theor. Phys. 112 (2004), 415. 
S. Mano and E. Takasugi, Prog. Theor. Phys. bf 97 (1997), 213. 
S. Mano, H. Suzuki and E. Takasugi, Prog. Theor. Phys. 95 (1996), 1079. 
[19] W. Schmidt, Class. Quant. Grav. 19 (2002), 2743. 
[20] L. Barack and A. Ori, Phys. Rev. D 64 (2001), 124003. 



